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ABSTRACT 

The formation of brown dwarfs (BDs) due to the fragmentation of early-stage 
proto-stehar disks undergoing pairwise encounters was investigated. High resolution 
allowed the use of realistic, flared initial disk models where both the vertical struc- 
ture and the local Jeans mass were resolved, preventing numerical enhancement or 
suppression of fragmentation. The results show that objects with masses ranging from 
giant planets to low mass stars can form during such encounters from disks that were 
stable before the encounter. The parameter space of initial spin-orbit orientations and 
the azimuthal angles for each disk was explored with a series of 48 simulations. For 
the types of interactions studied, an upper limit on the initial Toomre Q value of 2 
was found for fragmentation to occur. Depending on the initial configuration, shocks, 
tidal-tail structures and mass inflows were responsible for the condensation of disk gas 
into self-gravitating objects. In general, retrograde disks were more likely to fragment. 
It was also found that in fast encounters, where the interaction timescale was signifi- 
cantly shorter than the disks' dynamical timescales, the proto-stellar disks tended to 
be truncated without forming objects. 

The newly-formed objects had masses ranging from 0.9 to 127 Mj, with the ma- 
jority in the BD regime. These often resided in star-BD multiples and in some cases 
also formed hierarchical orbiting systems. Most of them had large angular momenta 
and highly flattened, disk-like shapes. With the inclusion of the appropriate physics 
these could reasonably be expected evolve into proto-BD disks with jets and subse- 
quent planet formation around the BD. The objects had radii ranging from 0.1 to 
10 AU at the time of formation. The disk gas was assumed to be locally isothermal, 
appropriate for the short cooling times in extended proto-stellar disks but not for 
condensed objects. An additional case with explicit cooling that reduced to zero for 
optically thick gas was simulated to test the extremes of cooling effectiveness and it 
was still possible to form objects in this case. Detailed radiative transfer (not studied 
here) is expected to lengthen the internal evolution timescale for these objects so that 
they spend considerable time as puffed-up, prolate ellipsoids, but not to alter the basic 
result that proto-BD disks can form during proto-stellar encounters. 



1 INTRODUCTION 



Observations suggest that brown dwarfs (BDs, 0.013 Mq < 
m < 0.075Mq) are very abundant in both the field and star 
clusters in our galaxy (|Chabrieij |2002| : iMartm et all l200ll : 
lLuhman|[2OO0l) . The origin of brown dwarfs, however, is still 
under debate. It has been suggested by s everal authors (e.g., 
IWhitworth et al.ll2007l : iBate et al.ll2002l ) that the final stel- 
lar mass may be closely related to the Jeans mass of the 
parent gas within the molecular cloud. Given that the typ- 
ical Jeans mass in a typical molecular cloud is of order one 
solar mass, this is reasonable for larger, hydrogen-burning 
stars. However, because BD masses are at least an order 
of magnitude smaller, it is difficult for BDs to form by di- 
rect collapse to pre-BD cores. Several scenarios have been 
proposed to circumvent this problem. 



One approach is to assume that the initial Jeans mass 
of the core is still a good estimate of the final mass from a 
direct collapse, but that denser or colder regions exist where 
the Jeans mass may appro ach the BD mass. For example, 
iPadoan fc Nordlundl (|2004l ) proposed that local density of in 
cores can be enha nced by turbulent flows ( "turbulent frag- 
mentation"), while lElmegreenl (|l999l ) suggested that certain 
types of molecular cooling can greatly decrease core temper- 
atures. In either case, formation of BDs follows the conven- 
tional star formation path. Hence, all properties associated 
with proto-stars, such as accretion disks and outflows, are 
expected. 

An alternat iv e p icture, first proposed by 
iReipurth fc Clarke! |2003), suggests that the because 
the Jeans mass continuously decreases during the collapse, 
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fragmentation can keep occurring repeatedly to produce 
smaller and denser objects until they are too dense to 
have efficient radiative cooling. The mass at this "opacity 
limit" is expected to be around mopacity ~ A/jupiter- Thus 
the opacity limit mass, rather than the initial core Jeans 
mass, would set the mass of the smallest fragments. In 
this scenario, it normally assumed that these small masses 
form first. These so-called embryos then accrete mass 
and ultimately form both sub-stellar and stellar objects 
in a scenario referred to as "competitive accretion". The 
accretion of the surrounding material onto embryos is 
assumed to be rapid so that embryos remaining in dense 
parts of the molec ular cloud would ultimately become 
stellar-mass objects. iBate et all l|2003l ) explored this sce- 
nario with a simulation of cluster formation from a 50 
Mq molecular cloud and found embryos could remain in 
the cloud undergoing rapid accretion or get ejected during 
violent interactions, resulting in a large number of low mass 
objects. More recent work indicates that a more realistic 
treatment of radiative energy losses can strongly suppress 
the formation of sub-stellar objects in such simulations 
ifPrice fc Bate 2009j) . 

The existence of extended disks around young stellar 
objects is supported by bot h observations and theory (e.g. 
lYorke fc Bodenheimeij Il999l ) . The disks are initially large 
and flared at the cl ass O-I stage (extending to several thou- 
sands AU (see, e.g.. iMundv et al.. 2000 ) and are observable 
at mm or sub-mm wavelengths. They evolve into smaller, 
thin proto-planetary disks at the class III stage that are 
detectable with infrared observations. Disk lifetimes are in- 
ferred to be about 10 Myr. Disks provide a natural environ- 
ment for forming sub-stellar objects because the Jeans mass 
in disks is much smaller than in cores. It may be possible to 
form brown dwarf mass companions through the fragmen- 
tation of margin ally stable (Q ~ 1), extended di sks around 
single stars (e.g.. IStamatellos fc Whitworthll2009h . Whether 
disks become cold and heavy enough to spontaneously frag- 
ment depends on the disk mass build up from the envelope 
and the processes that are driving accretion. Self-regulated, 
stable {Q > 1.5), extended disks are the outcome for a sig- 
nifica nt portion of the parameter space (|Vorobvov fc Basra 
I2OO9I ). 

The approach that will be investigated in this paper is 
encounter-triggered fragmentation of oth erwise stable pro to- 
stell ar disks. Numeri cal studies bv iWatkins et al.l (|l998al lbh 
and iLin et all (| 19981 ) have independently found that insta- 
bilities leading to bound companions can be triggered dur- 
ing encounters between stable proto-stellar d isks. Fragmen- 
tation was seen to occur within the disk (IWatkins et al.l 
Il998bl ). or in filamentary structures induced by the tidal 
force between the encountering disks l|Lin et al.l 19981 ). The 
resultant objects reported by both groups include sub-stellar 
mass objects. 

Since most stars form in clusters wit h a relatively 
small spread in ages (jGauvin fc Stromlll992l ), it is reason- 
ably likely that a proto-star will undergo a n encounter 
during the time that the d isks are large (e.g. iThies et al.l 
20051: IWatkins et al.1 Il998al ). Following the argument by 



eter at pericenter of 500 AU is about 30 percent. If sev- 
eral sub-stella r object s form during each encounter, as in 
IWatkins et al.l l|l998al ). a substantial numbers of BDs can 
be produced. In this work, we focus on the outcomes of spe- 
cific encounters rather than estimates of their likelihood. 

The prior simulation work failed to address several 
important issues. Firstly, the resolution used in the work 
cited above did not ensure that the Jeans mass was re- 
solved in the original disks. Thus it was not clear whether 
the objects might have condensed due to n umerical efi'ects 
(|Bate fc Burkertlll997l : lTruelove et al.lll997l ). 

Secondly, the range of circumstances in which gravi- 
tational instabilities can be triggered in a stable disk by an 
encounter were not explored. The Toomre criterion for insta- 
bility requires that that Q < Qait for thin disks to be unsta- 
ble to axisymmetric perturbations, whe re Q is the To omre 
parameter, locally defined according to l|Toomrelll964l ). 
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Thies et al.1 l|2005h . for a Plummer star cluster model with 
mass of 500Mq and a half-mass radius, Tio.s = 0.5 pc, the 
encounter probability in 5 x 10^ years (the approximate life- 
time for an early extended disk) with an impact param- 



where Cs is the sound speed, k is the epicyclic frequency 
(close to the angular speed Q for a Keplerian disk ) and S 
is the disk surface density. Analytical studies by iToomrd 
1 19641 ) and numerical studies of relatively thin proto- 
planetary disks (Boss 2 00ll : iMaver et aLllioM : iPickett et al.1 
l2003l : iGammij l200li ) have shown that Qcrit is in the range 
1-1.5 for various modes of fragmentation. If proto-stellar 
disks are stable according to this criterion then they tend 
to stay so indefinitely in the absence of mechanisms to re- 
distribute material. A non-negligible (vertical) thickness for 
the disks can be stabilizing, reducing the critical Q value 
for axisymmetric perturbations from 1 down to 0.6 — 0.7 
(|Kim et al.ll2002l ). The evolution of the Toomre parameter 
during dynamical interactions is a good indicator for the 
onset of instabilities and fragmentation and this has been 
explored here. It should also be noted that the initial disk 
models assumed in prior work did not resolve the disk ver- 
tical structure, which can destabilize the disk (by increas- 
ing the effective critical Toomre Q) and artificially enhance 
fragmentation. These issues motivated our use of detailed, 
vertically-resolved initial disk models. 

Lastly, there have been few predictions regarding the 
properties of the young BDs that could be produced during 
a proto-stellar disk encounter. Most observed young brown 
dwarfs show signatures of surrounding disks analogous to 
the disks around proto-stars, and some of these disks have 
been identified as accret ors through spectral line analysis 
(e.g. Ijavawardhana et al . 2003). The detection of BD accre- 
tion disks has been used as an evidenc e to support the tur- 
bulen t fragmentation scenario (see, e.g. lPadoan fc Nordlundl 
|2004 ) against dynamical formation mec hanisms, especially 
after the aforementioned simulation by iBate et aU (|2003l ) 
produced an practically no proto-BD disks. Ear l y enco unter 
simulations, s u ch as the ones by iLin et al.1 (Il998l ) and 
IWatkins et al] (Il998al lbh. also reported no proto-BD disks. 
However, these negative results are quite likely to reflect in- 
sufficient resolution or the impact of approximations such 
as replacing the fragments with a large "sink" particles 
(|Bate et al.l l2003l ). Studies with better resolved disks and 
fragments were thus required to address this question. 

The current study investigated encounter-induced BD 
formation using a series of high resolution hydrodynami- 
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cal simulations. This allowed the Jeans mass to be resolved 
until well after fragmentation initially took place and the 
vertical structure of the initial disks was also well resolved. 
The paper is organized as follows: in section [2] we describe 
the numerical method, the construction of the initial disk 
models and a study of resolution requirements. In section [3] 
we present the encounter parameter space being explored, 
including Q values, relative orientations of disks and varia- 
tions in encounter velocity. Section |3] describes several typ- 
ical simulations and discusses the conditions for encounter- 
induced fragmentation. Section [5] presents analysis of the 
properties of the resultant objects, including shapes, angu- 
lar momenta and orbital evolution. Section [7] summarizes 
the results and includes discussion of the observational im- 
plications and possible directions for future work. 



2.2 Radial Toomre Q Profile 

The choice of radial profile gave a nearly constant Q value 
between 100-500 AU, with a minimum value (Qmin) ranging 
from 1.6 to 2.1 depending on the case studied. The Q value 
rises in the inner 100 AU and the disk is very stable there. 
Test simulations of single disk evolution with different val- 
ues of Qmin were performed, and it was found that for our 
disk model values of Qmin below Qcrit ~ 0.8 were required 
for spontaneous fragmentation via spiral waves. This criti- 
cal value was lowe r than the ones see n in thin disk models 
(e.g., Qcrit = 1.4 in Mayer et al.l (j2003)), consistent with the 
results from lKim et al.l l|2002D that significant scale heights 
tend stabilize the disk. To rule out spontaneous fragmenta- 
tion as a factor during the dynamical interactions between 
disks, the Qmin values used were at least twice as the Qcrit 
for an isolated disk, i.e. Qmin ^ 1-6. 



2 METHODS 

2.1 Proto-Stellar Disk Model 

Observing the exact mass, density structure and kinematics 
of proto-stellar obj ects during the earliest stages is difficult 
(|White et al.ll2007 ). Using a hydrodynamical calcula t ion o f 
proto-stellar disk formation, lYorke fc Bodenheimeij (Il999l ) 
found that early disks contain a large fraction of the total 
mass of the system and collapse relatively slowly. Based on 
these results, we assumed that the proto-stellar disks were 
quasi-hydrostatic in order to construct a self-consistent disk 
density structure. 

We chose to simulate systems in the solar mass range 
with a single model for the disk-proto-star system where the 
disk mass was about 0.6 Mq, with a stellar mass of around 
0.5 Mq . The disk extended to about 1000 AU . This aspect of 
the set- up is similar to the disks simulated bv lWatkins et all 
l|l998al lbl). A power-law surface density profile E = Eq was 
chosen with p = 1.5 in 100-500 A U, consistent with obser- 
vations of early proto-stellar disks l|Andre et al.ll2000l ). The 
surface density peaked at 100 AU with T,max ~ 20.0 g/cmP' . 
Within 100 AU the gas was initially smoothly truncated. We 
relaxed our initial disks in isolation for about 10000 years, 
during which a small amount of material occupied the trun- 
cated region. Beyond 500 AU the surface density falls rapidly 
but smoothly, and at 1000 AU E ~ 0.03 g/cm?. The entire 
profile after relaxatio n is similar to the results from core col- 
lapse calculations bv lYorke fc Bodenheimerl l|l999l ). Most of 
the disk mass (0.5 Mq) is within 500 AU. 

It was assumed that the central proto-star was the 
major heating source for the disk, and that the disk re- 
radiates like a black -body, giving a di s k temperature pro- 
file of T(r) PC r~^/ ^ ('Mayer et al1l2004l : iPickett et al.ll2003l : 
IChiang &: Goldrdr h 1997) for the majority of the simula- 
tions. The profile smoothly transitions to an external, ambi- 
ent temperature of ~ 40 — 80 K. This profile is similar to the 
lYorke fc Bodenheimerl l| 19991 ) results. With the radial sur- 
face density and temperature profiles specified, the vertical 
structure was iterated to a self-consistent hydrostatic equi- 
librium taking into account both the gravity of the central 
star and the self-gravity of the disk. Accurate treatment of 
the gas self-gravity was important given the large disk mass 
relative to the star. 



2.3 Heating and Cooling Approximations 

A fixed temperature profile T{r) oc r'^^"^ for most of the 
disk was used in the majority of our simulations, which 
is often referred to as a locally isothermal equation of 
state (EOS). Heating sources other than the star, such as 
shocks, were not taken into account. For the current study, 
it was found that the isothermal EOS was a good approx- 
imation for the initial disk and also acceptable even after 
disk was heated by tidal structures and shocks (outside the 
clumps) because the disks are extended and initially opti- 
cally thin. Given the temperature of the gas, using the sur- 
face density profile described above and adopting the opac- 
ities from iD'Alcssio et al. (2001), the product of the disk 
cooling time tcooi and angular rotation rate Q. was calcu- 
lated and Q. tcooi = 0.3 — 0.7 beyond 50 AU for the ini- 
tial disk, which is comfortably within the cooling require- 
ment for the dis k to fragment Q. tcooi < 12 (H ice et al.ll2005l : 
Ijohnson fc Gammie 2003| ) before the disk is dense enough 
for the isothermal assumption breaks down. However, after 
the disk fragments and the resultant objects became grav- 
itationally bound, the final densities of these objects were 
typically orders of magnitude larger than the initial values 
and the objects were thus optically thick. Thus our locally 
isothermal approximation was initially good but later be- 
came an upper bound on how efficient cooling should be. 
To further examine the importance of cooling in fragmenta- 
tion, a case was also investigated with explicit cooling that 
became negligible in the optically thick regime. This second 
assumption is a lower bound for the cooling in a real disk 
as it disallows even the slow leakage of radiation from dense 
gas. The results of this case are described in section [S] 



2.4 Simulations 



We u sed the TreeSPH code GASOLINE (|Wadslev et all 
l2004h to model the disks during encounters. Around 200,000 
particles were used for each 0.6Mq disk, with each particle 
representing ~ 1 Earth mass. This resolution was chosen 
to satisfy the condition that the lo cal Jeans mass should 
be resolved by at least 2Nneighbor iBate fc Burkertlll997l ) 
to prevent numerical fragmentation, where Nneighbor is the 
number of neighbor particles in SPH (Smoothed Particle Hy- 
drodynamics) simulation. The number of particles ensured 
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that the Jeans mass was not-only well resolved in the initial 
condition, but also well-resolved when density was enhanced 
during an encounter (when the Jeans mass decreased), up 
to a density p ~ 10"^" g cm~^ (about 10^ times the ini- 
tial maximum value in the disk). Resolving the Jeans mass 
within the fragments, however, was not attempted due to 
high computing costs (A'^ > 10^ particles usually required). 
Therefore, our choice of particle numbers ensures physically 
correct modeling of the initial fragmentation in the proto- 
stellar disk, but cannot ensure that any secondary fragmen- 
tation of the clumps themselves are physical. 

The initial disk model also resolved the scale height at 
the mid-plane very well at 100 AU (where the local smooth- 
ing length was about 1.9 AU and the scale height was about 
6 AU). Within 50 AU the scale height of the disk was less 
well resolved and the orbital time was shortest so that artifi- 
cial viscous evolution could have been a concern. The artifi- 
cial viscosity was moderated using the Balsara switch which 
helps to suppress unwanted shear effects. The measured vis- 
cous time scale in the code was about 10* to 10^^ years, much 
longer than the dynamical scale of the encounter, t ~ 10^ 
years. Hence, the effect due to the viscous evolution was neg- 
ligible. Appendix |X] presents as exploration of the effects of 
unresolved Jeans mass and unresolved disk scale height. 

The gravitational softening, e, for gas particles was 0.2 
AU, one tenth of the smallest particle spacing in the initial 
condition. Thus dense regions up to 1000 times the initial 
densest part were well resolved. Since a larger e usually sup- 
presses gravitation forces, the choice of gravitational soft- 
ening will tend to inhibit fragmentation once the particle 
separations are comparable to or smaller than the soften- 
ing. In all tests, including those done for this work and in 
other work to date, there have been no indications that the 
Gasoline code suffers from artificial fragmentation as a con- 
sequence of the gravitational softening being smaller than 
the local SPH smoothing length. The fixed softenings en- 
sure accurate gravitational dynamics at the cost of having 
somewhat conservative time steps in lower density regions. 

The central proto-star was represented by a single 0.5 
solar mass sink particle for which the softening was set to 1 
AU. The sink acted within a 1 AU radius around the star 
to consume gas falling into the innermost regions and thus 
avoid the small time steps required for gas orbiting very 
close to the star. However, no objects formed during the 
simulation were replaced with sink particles, ensuring that 
any disks surrounding the newly formed, sub-stellar object 
were modeled directly. 

3 PARAMETER SPACE 

3.1 Locally Isothermal Simulations 

In each pairwise disk encounter, the two disks were identical 
in size, mass and Toomre parameter with an initial impact 
parameter of 1000 AU. Due to gravitational focusing during 
relatively slow encounters the minimum distance between 
two stars could become substantially less than the initial 
impact parameter. The initial disk separation was 6000 AU. 
This system was initially bound if the relative velocity of 
the encounter is less than 0.8 km/s and the final impact 
parameter is a function of the relative velocity. This study 
explored six dimensions in the overall parameter space: 



(i) The initial Toomre parameter profile characterized us- 
ing the minimum value, Qmin; 

(ii) The geometric configuration of the encounter, includ- 
ing the angles between the two disk planes and the angles 
between the disk planes and the plane of the orbit of the 
two disks (a total of four independent parameters); 

(iii) The relative velocity of the encounter. 

Mass, disk size and impact parameter variations were 
not included to limit the study to a manageable size. Smaller 
impact parameters are possible but considerably less likely. 
These variations may be explored in future work. 

A list of all the simulations is shown in Table [T] The 
first 3 cases (Ql, Q2, Q3) were carried out to study the effect 
of the Toomre parameter on encounter-triggered fragmenta- 
tion with different initial minimum Q values of 1.6, 1.8 and 
2.1, respectively. A higher Q profile means a more stable 
initial condition, for which more effort is required to trig- 
ger fragmentation. We selected an encounter configuration 
where fragmentation was likely (coplanar, retrograde disks) 
to first explore the threshold in initial Toomre Qmin, above 
which dynamical interactions would become insufficient to 
trigger fragmentation. This study helped to constrain the 
region of interest with in the par ameter space. 

Previous work bv lLin et al.l |l998) indicated that frag- 
mentation can ta ke place in tida l structures forming during 
encounters while IWatkins et al.l l|l998al lbl) suggested disks 
would fragment. The degree of tidal distortion and angu- 
lar momentum transport has been shown to depend on the 
initial disk spin orie ntations (with respect to the orbital an- 
gular momentum) (|Toomre fc Toomrd Il97 j ). Various disk 
spin orientations were used in our simulations including pro- 
grade cases (where the angle between the angular momen- 
tum vector of the disks and that of the orbit was less than 
90 degrees), retrograde cases (where the spin-orbit angle was 
greater than 90 degrees) and combinations of the two. 

In most simulations the disks were not close to the same 
plane, as would be expected for randomly oriented disks. 
However, coplanar encounters (where the spin-orbit angles 
are zero) were also simulated as special cases to investigate 
the shocks generated in such encounters and the possibility 
of shock fragmentation. In every case the orbital plane was 
the x-y plane so that the orbital angular momentum was in 
the direction of the z-axis. 

In a completely generic encounter, as shown in Figure [TJ 
the disks can be oriented randomly. The configuration of 
each disk in the encounter can then be characterized by two 
angles: the angle between the spin angular momentum of 
each disk and the positive z direction, 6, and the angle be- 
tween the positive x direction and the projection of disk 
angular momentum onto the x-y plane, <j), which also deter- 
mines where that disk intersects the x-y plane. The initial 
disk velocities are always along the x-axis and the initial 
impact parameter is due to an offset in the y-direction. 

Cases Confl-48 explored the {9i, (f>i, 62, 4>2) parameter 
space. For a single disk, the angle parameters were selected 
such that the 9 space (from to 180 degrees) was divided 
evenly into three ranges and the (j) space (from to 360 
degrees) into four for a total of 12 cases. Naively there would 
have been 144 possible combinations of disks but due to the 
symmetry of the encounters and the equivalence when the 
two disks were exchanged in space, there were 42 unique 
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Figure 1. A generic encounter configuration between two proto-stellar disks. Tlie centres of mass orbit in tlie x-y plane witli an initial 
velocity solely in x-direction. The angular momentum of an individual disk has a vector direction that differs from the orbital angular 
momentum (the z-axis) by an angle 9. The disk angular momentum vector projected onto the x-y plane makes an angle (j) with the x-axis. 



collisions to simulate. An additional six cases were simulated 
where at least one disk was exactly in the orbital plane, to 
investigate the effects of coplanar encounters. 

The cases with labels starting with "Vel" refer to two 
series of runs used to examine the constraints due to the ini- 
tial relative velocity on disk fragmentation for prograde and 
retrograde encounters. Relative encounter velocities in the 
range of 0.8 kms~^ to 6.0 kms~^ were used in the simula- 
tions to investigate the effects of relative velocity on disk 
fragmentation. The inferred disk velocities are consistent 
with the typical velocity disp ersion in stellar clusters (a few 
kms~^, See e.g.. lBate et al.li2003i . and references therein) 



4 NUMERICAL RESULTS OF THE 
PARAMETER SPACE STUDY 

4.1 Varying Toomre Q parameters 

To explore the impact of the initial minimum Toomre pa- 
rameter we investigated 3 cases, all of which have {9, 0) = 
(180,0) and relative v — 0.8 km/s. These systems are only 
marginally unbound. These collision parameters were cho- 
sen to ensure a strong interaction (detailed discussion in 
Section 14. 2p so that if no fragments form in these calcula- 
tions then they are also not likely under other configurations 
and encounter speeds. The two disk are coplanar since they 
have the same 6 and <j>. 

Case "Ql" (Qmin=1.6) is depicted in Figure [21 When 
the physical impact began, the gas was shock compressed 
and quickly formed a "shock layer" between the two disks 
where the density was substantially enhanced. At 24,000 
years after the encounter, the shock fragmented to form a 
7 Mj clump at about 223 AU away from the nearest star 
(label "a" in Panel 5, Panel 7). Also, due to its retrograde 
configuration, the angular momentum of the disk spin and 
the orbit are opposed and thus the rotation speed of the 
gas decreases. As a result, a large amount of gas spirals into 
the inner disk (within ~ 100 AU in this case), while gas in 



the outer disk dissipates. The azimuthally averaged surface 
density profile of the lower-left disk in Panel 4 of Figure [2] is 
shown at different times in Figure [3] and the corresponding 
evolution of azimuthally averaged Q profile is shown in Fig- 
ure|4] At t = 23,000 years, the surface density had increased 
significantly, and at r ~ 10 — 20 AU, it had almost doubled 
its value from earlier times (t — 6,000 years, well before the 
physical impact of gas materials). However, Q was still very 
high in the inner 100 AU because of the fixed high temper- 
ature in that region (Figure |4|. Thus it was the presence 
of the shock layer at around 200 AU that locally increased 
the surface density and lowered the Toomre Q down to unity 
(Qmin ~ 1, Figure ID) allowing a strong spiral feature to form 
which is gravitationally unstable. At t = 24,000 years, a pair 
of close clumps formed at r ~ 130 AU with masses of about 
12 and 14 M,j (labeled "b" in Panel 5, Panel 8), and at 
about 150 AU another gas blob started condensing (a mod- 
erate extremum near the 130 AU clump in the plot of the Q 
profiles, labeled "c" in Panel 5), which eventually became a 
bound object (labeled "d" in Panel 6). Gas inflow increased 
after formation of these objects. At t = 24,400 years, frag- 
mentation within 20 AU was finally triggered in one of the 
disks and 4 brown-dwarf mass clumps were formed (boxed 
region in Panel 6, Panel 9). However, these objects were not 
included in our statistics (Section [S} because the density 
within 20 AU at this time step exceeded W~^'^g cm~^ and 
became optically thick (cf. Section O, so that an isothermal 
cooling assumption was no longer appropriate. The simula- 
tion was stopped at t = 24,400 years because the time step 
reduced dramatically after the first condensation formed. 
Analysis based on the energy and momentum at the end of 
simulation indicates that amongst the 4 objects formed out- 
side 100 AU, one is unbound and likely to leave the system, 
one has a large semi-major axis (a = 879.3 AU) and the 
other two have orbits with semi-major axes less than 100 
AU. 

The encounter "Q2" has same parameters as "Ql" ex- 
cept the initial Qmin ~ 1-8. The disk underwent similar 
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Table 1. Simulation Parameters 



Case 


Toomre Qniin 


(t'l><Pl) (deg) 


(('2,92) (deg) 


V (km/s) 


Nl 


AT 

1\2 


Ql (VeLretrol) 


1.6 


(ISO, 0) 


(180, 0) 


0.8 


3 (1) 


(4) 


Q2 


l.S 


(ISO, 0) 


(180,0) 


0.8 





(1) 


Q3 


2.1 


(180,0) 


(180,0) 


0.8 








Confl 


1.6 


(180,0) 


(180,0) 


2.0 


(1) 


(1) 


Conf2 


1.6 


(0,0) 


(0,0) 


2.0 


3 


4 


ConfS {Vel.pro2) 


1.6 


(0,0) 


(45,0) 


2.0 


1 


2 


Conf4 


1.6 


(180,0) 


(135,0) 


2.0 


1 


(1) 


Conf5 


1.6 


(180,0) 


(45, 0) 


2.0 


2 (1) 





Conf6 


1.6 


(150, 110) 


(150, 110) 


2.0 


1 (1) 


2 (1) 


Conf7 


1.6 


(150, 110) 


(150, 250) 


2.0 





4 


ConfS 


1.6 


(150, 110) 


(150, 290) 


2.0 





3 (1) 


Confg 


1.6 


(150, 70) 


(150, 70) 


2.0 


2 


3 


ConflO 


1.6 


(150, 70) 


(150, 110) 


2.0 





3 (1) 


Confl 1 


1.6 


(150,70) 


(150,250) 


2.0 


4 


3 (1) 


Confl2 


1.6 


(150,70) 


(150,290) 


2.0 


1 (1) 


2 (1) 


ConflS 


1.6 


(30, 110) 


(150, 110) 


2.0 





9 (1) 


Confl4 


1.6 


(30, 110) 


(150,250) 


2.0 


1 (1) 


(1) 


Confl5 


1.6 


(30, 110) 


(150, 290) 


2.0 





5 (1) 


Confie 


1.6 


(30, 70) 


(150, 70) 


2.0 





3 (2) 


Confiy 


1.6 


(30, 70) 


(150, 110) 


2.0 





6 (1) 


ConflS 


1.6 


(30, 70) 


(150, 250) 


2.0 








ConflQ 


1.6 


(30, 70) 


(150,290) 


2.0 





3 (2) 


Conf20 


1.6 


(30, 110) 


(30, 110) 


2.0 


3 (1) 





Conf21 


1.6 


(30, 110) 


(30,250) 


2.0 








Conf22 


1.6 


(30, 110) 


(30,290) 


2.0 


3 


2 


Conf23 


1.6 


(30, 70) 


(30, 70) 


2.0 


1 (1) 


1 (1) 


Conf24 


1.6 


(30, 70) 


(30, 110) 


2.0 








Conf25 


1.6 


(30, 70) 


(30, 250) 


2.0 





3 


Conf26 


1.6 


(30, 70) 


(30, 290) 


2.0 


1 





Conf27 


1.6 


(30, 110) 


(90, 135) 


2.0 





2 


Conf2S 


1.6 


(30, 110) 


(90,225) 


2.0 





1 (1) 


Conf29 


1.6 


(30, 110) 


(90,315) 


2.0 








Conf30 


1.6 


(30, 70) 


(90,45) 


2.0 





3 (1) 


ConfSl 


1.6 


(30, 70) 


(90, 135) 


2.0 





5 


Conf32 


1.6 


(30, 70) 


(90,225) 


2.0 





1 


Conf33 


1.6 


(30, 70) 


(90,315) 


2.0 





3 


Conf34 


1.6 


(90, 135) 


(150, 110) 


2.0 





4 (5) 


Conf35 


1.6 


(90, 135) 


(150, 250) 


2.0 


5 (1) 


5 


Conf36 


1.6 


(90, 135) 


(150, 290) 


2.0 


3 (3) 





Conf37 


1.6 


(90, 45) 


(150, 70) 


2.0 





2 (1) 


ConfSS 


1.6 


(90, 45) 


(150, 110) 


2.0 


2 


3 


Conf39 


1.6 


(90,45) 


(150, 250) 


2.0 


4 (1) 





Conf40 


1.6 


(90,45) 


(150, 290) 


2.0 





5 


Conf41 


1.6 


(90, 135) 


(90, 135) 


2.0 








Conf42 


1.6 


(90, 135) 


(90,225) 


2.0 


2 (1) 





Conf43 


1.6 


(90, 135) 


(90,315) 


2.0 








Conf44 


1.6 


(90,45) 


(90,45) 


2.0 








Conf45 


1.6 


(90, 45) 


(90, 135) 


2.0 








Conf46 


1.6 


(90, 45) 


(90, 225) 


2.0 








Conf47 


1.6 


(90, 45) 


(90,315) 


2.0 


1 


4 


Conf4S 


1.6 


(0,0) 


(90, 0) 


2.0 





3(1) 


VcLprol 


1.6 


(0,0) 


(45, 0) 


0.8 


0(1) 





Vel_pro2(Conf3) 


1.6 


(0,0) 


(45, 0) 


2.0 


1 


2 


Vel_pro3 


1.6 


(0,0) 


(45, 0) 


4.0 








VeLpro4 


1.6 


(0,0) 


(45, 0) 


6.0 








Vel_retrol(Ql) 


1.6 


(180,0) 


(180,0) 


0.8 


3 (1) 


0(4) 


Vel_retro2 


1.6 


(ISO, 0) 


(180,0) 


2.0 


0(1) 


0(1) 


Vel_retro3 


1.6 


(ISO, 0) 


(180, 0) 


4.0 


1(1) 





Vel_retro4 


1.6 


(ISO, 0) 


(180, 0) 


6.0 









The pairs (6i, t^i) and {62, 4>2) denote the angles between the disk angular momentum vector 
and the angular momentum vector of the orbit and the angular momentum vector and the 
vector direction of the initial velocity difference in the orbital plane for disk 1 and disk 2 
respectively, v denotes the relative encounter velocity. The number of objects formed for 
each disk are shown in the last two columns. If objects were formed within 50 AU then the 
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Figure 2. A coplanar encounter between two retrograde disks witli Qmin = 1-6 . Panels 1 to 6 are snapshots of tlie simulation, with 
the time indicated at the left top of each panel. The disk were rotating clockwise in all of the panels. The color bar under each panel 
gives the logarithm (base 10) of the density in units of Mq / AU^ . In panels 4 and 5 only the inner dense regions were included and a 
different density range was used as indicated. Panels 7 and 8 are zoomed-in pictures of objects formed from a shock layer and due to disk 
fragmentation (a and b in panel 5), respectively. Panel 9 is the zoomed-in picture of the boxed region in panel 6, showing fragmentation 
of the inner disk. 
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100 
r(AU) 



Figure 3. Evolution of tlie surface density profile for the case in 
Figure [2] Solid line: Surface density profile of one of the encoun- 
tering disks at 6,000 years. Dashed line: Surface density profile 
for the same disk at 23,000 years. Dotted line: Surface density 
profile at 24,000 years. The sharp peak at 130 AU is produced by 
two close objects forming due to disk fragmentation, the one at 
225 AU originates from shock layer fragmentation. The moderate 
peak in between indicates ongoing gas condensation. 




Figure 4. The Toomre Q profile for the same disk as in Figure 
[3] at 6,000, 23,000 and 24,000 years (solid, dashed, dotted Une, 
respectively). The dot-dashed line indicated Q = 1.0 



physical processes such as formation of a layer of shock- 
compressed gas, inflow of gas and disk truncation after the 
encounter. Nevertheless, the disks were more stable to the 
perturbation and only one of the disks developed an m = 2 
mode gravitational instability at a relatively late time step 
(t = 25,300 years), which then fragmented and formed only 
one bound object at radius of 30 AU. The other disk was 
truncated but remained stable. The gas in the shock layer 
did not condense to form an object. Instead, it dispersed af- 
ter the disks moved apart. The lowest Q value at t = 24600 
years (the time step just before fragmentation) for the un- 
stable disk was 1.1 and for the stable disk was 1.7, close to 
the original minimum. 

The last case in this series "Q3" used initial disks with 
Qmin = 2.1. The encounter disrupted both disks without 
forming any objects. The two stars happened to capture each 



other and formed a close binary system with a separation of 
about 20 AU. Most of gas was dispersed, with the bound gas 
at the end of simulation forming a circumbinary disk with 
a radius of 250 AU at 31000 years. The circumbinary disk 
had a much lower surface density compared to the initial 
disks. Assuming the gas was still in Keplerian rotation, the 
minimum Toomre Q was very high (Qmin = 6.2). The disk 
was thus very stable against spontaneous fragmentation and 
no clumps were expected to form in future evolution. 

The cases described here showed that instability and 
fragmentation can be triggered in an initially stable disk by 
pairwise encounters. However, the initial minimum Toomre 
parameter can not be higher than 2.1 for instabilities to oc- 
cur. From the evolution of the Q profile in Figure|4]it can be 
seen that lowering Q down to unity is crucial for fragmen- 
tation. Two competing processes, gas inflow and dispersion, 
proceed at the same time in an encountering disk. In a high 
Q disk the gas is hotter and an encounter is more likely 
to accelerate the gas dispersion rather than triggering gas 
inflow as in the relatively low Q (but still initially stable) 
cases. Indications are that early stage proto-stellar disks are 
usually large, heavy and relatively cold, with temperatures 
not far from those of the molecular cloud core (~ lOK). So 
gas inflow can be the key mechanism in encounter-induced 
disk fragmentation. 

Note that the local enhancement of density in the 
shocked gas can lower the Q parameter significantly and 
induce disk instability before the inflowing gas reduces Q 
in the inner, hotter region, making the fragments form rel- 
atively far from the stars. This is especially clear for copla- 
nar encounters discussed in this section. In more likely, non- 
coplanar cases (cf. Section |4]2} the shocked gas usually does 
not form a well-deflned layer between the disks, but it re- 
mains an important factor driving disk instability. In the 
case "Ql" the shock layer itself also fragmented directly. 
The object formed this way was not associated with spiral 
structure in the unstable disk (panel 5 in Figure [2} . Frag- 
mentation of shock layers rarely occurs and was restricted to 
coplanar, low-velocity, low-Q cases. With higher disk tem- 
perature (as in "Q2" ) or higher velocity the shock layer usu- 
ally dissipates when the two stars pass periastron without 
forming objects. 

4.2 Varying the Encounter Configurations 

The parameter space of encounter conflgurations was ex- 
plored with 48 encounters covering the (9, (f>) parameter 
space. The effect of the relation between disks' spin and or- 
bital angular momentum in z direction can be see by varying 
the disks' 9 parameters. This variation qualitatively changes 
the interactions, including the mechanisms of fragmentation. 

There is a large stochastic element to the evolution of 
the disks and the formation of clumps. In this sense, the pre- 
cise number of objects formed in a single simulation should 
not be examined too closely. For example, a small perturba- 
tion in the initial particle placements that is within the noise 
in the initial surface density profile is sufficient to change the 
outcome in terms of the precise number and placements of 
clumps. For example, several cases are symmetric from the 
point of view of the two disks but have non-symmetric out- 
comes. Symmetry in the conflgurations is revealed by the 
mappings x,y,z — x, —y, z which implies Oi = 62 and 
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i^i = 4>2 (such as for the coplanar cases and "Conf6", for 
example) and x,y,z — ^ ~x, —y, z which gives d\ — 62 and 
4>i = 4>2 + 180 (e.g. "Confll"). These initial states were 
created using absolute offsets {x ^ x + /S.x) of two iden- 
tical disks that have small non-axisymmetric perturbations 
due to the glass initial particle state. Thus the small pertur- 
bations in the particle configurations are not symmetric in 
these cases. These difference magnify during the course of 
the encounter and ultimately result in different numbers of 
objects in each disk as can be seen in the object numbers for 
each disk in symmetric cases in table [T] For example disk f 
in "Conf6" made 2 objects while Disk 2 made 3. The vari- 
ations are consistent with a Poisson-type random process. 
In most of our analyses we look at outcomes averaged over 
several encounters. 




100 

Angle 8 (Degree) 



4..2.1 Formation of Tidal Structures in Prograde Disks 
and Fragmentation 

Encounters between retrograde disks tended to reduce the 
disk spins and cause gas infiows, increasing the surface den- 
sity of the inner disks and reducing the Toomre Q parameter 
(Figure [2]). Prograde encounters are more complex as seen 
in Figure [5] showing snapshots of encounter "ConfS" . The 
two disks have {9,(f)) = (0,0) and (45,0), respectively. The 
encounter had a relative speed 2.0 km/s and the initial min- 
imum Q of the disks was 1.6. Gas between the two disks 
was shock-compressed (Panel 2 in Figure [5l , though it did 
not form a layer or fragment. In prograde cases spin-orbit 
resonances can occur. Gas in the resonance region is contin- 
uously perturbed promoting angular momentum transport 
that simultaneously enhances the density of the inner disk 
while the outward gas flows collect a substantial portion of 
the disk mass into filamentary tidal structures (Figure [S] 
Panel 4), analogous to the tidal tail structures generated 
during the interaction between two prograde disk galaxies. 
In the case shown, at 14,000 years gas began condensing 
from a tidal structure of one disk (on the right in Panel 
4) at about 320 AU from the star (labeled "a" in Panel 
4), which later formed a bound object of mass ~ 20 Mj. 
A tidal structure was also present in the other disk, where 
there was a 45 degree angle between the spin and orbital 
angular momentum vectors. At 15,000 years, it condensed 
into two objects with masses of 9 M,j and 70 M,j at 60 and 
110 AU from the star, respectively (labeled "b" in Panel 5 
Figure [SJ. The simulation was terminated at 16,000 years. 
Analysis based on the energy and momentum of the objects 
found that all the 3 objects remained bound to the system 
but their orbits were large: 532 AU for the one that con- 
densed first, and 90 and 307 AU for the two condensed in 
the later time step. 

Fragmentation of tidal structures was seen in several 
prograde disk encounters, especially when the disk was 
closely ali gned with the or bital plane. This result is consis- 
tent with iLin et al.l (| 19981 ). in which a disk-disk encounter 
triggered formation of a large tidal tail and an object con- 
densed from it. Objects formed this way tend to have large 
orbital radii. For example, in the case"Conf2", tidal struc- 
tures from both disks condensed and formed 7 objects in 
total by the end of simulation, all of which were unbound to 
both stars. 



Figure 6. The number of objects formed in a disk per encounter 
as a function of 6. Each box corresponds to a single simulated disk 
that produced the number of objects indicated on the y-axis. The 
boxes have been spaced apart in the x-direction for clarity. The 
average and variance are indicated with solid and dashed lines 
respectively. 



4-2.2 Prograde vs. Retrograde Encounters 

Tidal structures due to spin-orbit resonances are not as ef- 
ficient at producing objects as the inner disk fragmentation 
seen in retrograde disks. On the other hand. IWatkins et al.l 
(1998a) has argued that star formation triggered by col- 
lisions of molecular clouds tends to result in prograde- 
prograde coplanar encounters so they may be more common. 

In Figure |B] we plot the number of object formed in a 
disk per encounter as a function of 9, the angle between the 
individual disk spin vector and the positive z direction. We 
chose the homogeneous sample from "Conf6" to "Conf48" in 
Table [1] excluding the unlikely coplanar encounters, (p — 
encounters and encounters with disks residing in the orbit 
plane. The cases chosen all have same initial Toomre param- 
eter and relative encounter speed. 

The average number of bound and unbound objects 
formed within that disk increases with increasing 6 for that 
disk. This is the case regardless of the configuration of the 
other disk in the interaction. Gas infiow efficiently adds mass 
to inner disks which fragment whereas in prograde disks, the 
mass transfer usually send substantial amounts of gas out- 
ward rather than inward. The condensation of tidal struc- 
tures did not occur frequently, because in most of the cases 
they usually dissipated before anything could form. 



4.2.3 Plane-Parallel Encounters 

No clear correlation was found between the productivity of 
a disk and the disk's own (p parameter. However, no frag- 
mentation would happen if the two disk planes are parallel 
(or nearly parallel) initially but not coplanar. 60% of the 
encounters in the "Conf" series which do not produce any 
clumps fall into this category. The disks were simply trun- 
cated in these cases. The final size of the disk is about 300 
to 500 AU depending on the initial configuration. 
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t= 10,000 years 








Figure 5. Snapshots of the prograde-prograde disk encounter "ConfS". Time is indicated on the left top of each panel. The disks are 
counter-clockwise rotating in this view. The color bar under each panel gives the logarithm (based 10) of the density in units of Mq /AU^ . 
Different density ranges are used to clearly show the clumps. 



4.3 Varying Encounter Velocities 

We investigated the effect of varying the relative velocity 
with two series of simulations, one containing non-coplanar 
encounters with both disks prograde, and one with coplanar 
encounters with two retrograde disks. In each series, the rel- 
ative velocity of the encounters is repeated along with the 
number of object formed in Tabled 

With the exception of a special single case ("VeLprol", 
see below), the number of clumps formed in an encounter 
generally decreases with increasing velocity. No clumps 
formed in prograde disks when the relative encounter ve- 
locity was above 4.0 km/s. This result can be qualitatively 
understood by comparing the interaction timescale with the 
dynamical timescale of the disks. Though gravitational in- 
teractions occur regardless whether or not there is a physical 
impact, these simulations indicate that significant changes 
in the density profile occur only when the disks are close 
enough to physically impact. The interaction timescale is 



Table 2. Number of objects formed in the simulations with dif- 
ferent encountering velocities 



Case 


Velocity 


Number 


Fragmentation 




(km/s) 


of Clumps 


Mechanism 


Vbl_prol 


0.8 


1 


Circumbinary disk 


Vel_pro2 


2.0 


3 


Tidal structure and disk 


VeLproS 


4.0 







Vel_pro4 


6.0 







Vel_retrol 


0.8 


8 


Shock layer and disk 


Vel_retro2 


2.0 


2 


Disk only 


Vel_retro3 


4.0 


3 


Disk only 


Vel_retro4 


6.0 


1 





approximately. 



where d^s is the twice the disk diameter (4000 AU in this 
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case), which gives the rough separation for physical impact. 
V is the relative velocity of the encounter. The dynamical 
timescale within each disk is a function of radius r, given 
by, 



(3) 



where G is the gravitational constant and M* is the mass 
of the star. The encounter timescales tint were 2.4 x 10*, 
9.5 X 10^, 4.7 X 10^ and 3.1 x 10^ years for relative speeds 0.8, 
2.0, 4.0 and 6.0 km/s, respectively. The dynamical timescale 
tdyn was about 1.1 X 10* years at 400 AU and 3.9 x 10^ year 
at 200 AU. The implication is that there was insufficient 
time during fast encounters, where tint < ^dyn to apply the 
torque that transfers angular momentum to change the disk 
profiles and causes fragmentation. This is more restrictive 
for prograde disks where tidal structures form further out 
at distances similar to 400 AU. 



M/Mj 



Figure 7. The mass distribution of 146 objects formed outside 
50 AU in the 48 "Conf" simulations. The bin size is 5 A/j. 



4.4 Captures and the Formation of Binary Stellar 
Systems 

In the case "VeLprol", star capture occurred and the two 
encountering stars formed a binary system. Both disks were 
prograde and large tidal tail structures formed when the 
disks passed periastron. The tidal structure dissipated over 
time and carried away the excess angular momentum from 
the system, enhancing the formation of a close binary. The 
separation of the two components was about 90 AU, signif- 
icantly smaller than the periastron of the trajectory (about 
300 AU taking into account the effect of gravitational focus- 
ing). The circumbinary disk was about 500 AU in radius. 
At a later time it also fragmented to form a 20 Jupiter mass 
object. 

Another case where stellar capture occurred is the case 
"Q3" (as described in Section|4]l}. Since the starting disk Q 
parameter is 2.1, the circumbinary disk became hotter and 
less massive compared to the standard cases. No fragments 
formed at the end of the simulation. The initial disks are 
retrograde and the encounter did not induce tidal structures. 

The two cases above are the only ones where stellar cap- 
ture occurred. Note that both of them have rather conser- 
vative choices of encounter velocities (0.8 km/s), with which 
the system are marginally bound at the initial state. 



5 PROPERTIES OF THE OBJECTS FORMED 

Most of the encounters produced clumps with the excep- 
tion of those with high relative velocities or very stable ini- 
tial disks (e.g. Qmin > 2.1). As described in the last sec- 
tion, bound objects can form from fragmentation of shock 
layers, instabilities within a disk and in tidal structures. 
Since the timescale of disk encounters (about 10,000 years) 
is short and encounters may happen when there is still a 
lot of material obscuring the stars, encounters in progress 
are unlikely to be observed whereas the sub-stellar objects 
formed may persist indefinitely. Though these objects may 
evolve substantially from their time of formation, particu- 
larly in terms of luminosity, simulated clump properties such 
as mass, angular momentum and orbital kinematics can be 
directly compared with properties of observed objects. In 



the initial state, the disks were well resolved outside 50 AU. 
However, as the simulations progressed, additional matter 
accumulated in the inner disks improving the resolution el- 
ements per scale height. On the other hand, objects with 
50 AU are likely to experience further evolution such as ac- 
cretion and migration into the star. In addition, work by 
iRafikovl |2007) has suggested that disks this close to a star 
are optically thick and unlikely to cool fast enough to frag- 
ment. For these reasons, we have excluded objects within 
50 AU from the sample discussed below. However, the ob- 
jects formed within 50 AU are not statistically different and 
properties such as the relative mass distribution would not 
be significantly altered if these objects were included. A to- 
tal of 191 objects formed in our simulations (excluding the 
high Q cases) with 144 outside 50 AU. 



5.1 Mass Distribution 

The object masses range from 0.9 Jupiter masses (Mj) to 
127 Mj, extending into the stellar regime. The mass distri- 
bution is plotted in Figure [T] It was found that the number 
of objects generally decreases with increasing mass, con- 
sistent with the observed sub-stellar initial mass function 
(IMF). The simulated population does not include objects 
much below 10 Mj). In this regime there are 1000 or fewer 
particles per object. The numerical gravitational softening 
(0.2 AU) could also inhibit the formation of smaller objects. 
Thus these simulations are unable to probe whether there is 
a physical mass cut-off in this range. 

Observationally, the sub-stellar IMF in open clusters 
can be well described by a power-law dN/dM oc M~". 
For young open clusters such as Pleiades, A Orionis and 
cr Orionis, the exponent a is close t o 0.6 at low m asses 
with an uncertainty of order 0.1 (Moraux et al.l 120031 : 
iBarrado v Navascues et al.|[200i : ICabaUerq et al.ll2007l ') . An 
universal IMF was proposed by I KroupaT 2001 ) in which 
a = 0.3 ± 0.7 in the range O.OIMq < M, < O.OSMq. In 
this work, we fit the mass distribution of the sub-stellar re- 
gion (M < O.OSMq) and the best fit gives a = 0.6 with 
uncertainty 0.15 (Figure[S)l, consistent with the observed 
values. Above the stellar boundary the number of objects 
declines more rapidly with mass and a break is apparent 
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Figure 8. The logarithmic mass distribution of the objects. The 
error bars were calculated assuming Poisson errors oc V Afbin • 
Stars: sub-stellar objects (M ^ 80Mj); Diamonds: low mass stars 
{M > SOMj); Dashed line: linear fit of the logarithmic distribu- 
tion of sub-stellar objects dN/dM oc M~°', with a = 0.6 

near 0.08 Mq (or about 0.1 of the central mass) which i s 
also consistent with the trend reported by (lKroupall200ll ). 
There are only 9 objects in total that are beyond 0.08 Mq 
and thus the statistics are poor. More simulations, includ- 
ing those with varying central star masses, are required to 
explore the resulting mass function in the low mass, stellar 
regime. 

5.2 Shapes and Angular Momenta 

Most of resultant clumps in our simulation are not spherical 
but have highly flattened, disky shapes (height to radius ra- 
tio ~ 0.1), as shown in Figure|51 The radii of the disks range 
from 0.1 AU to 10.0 AU. The lower bound is likely to be as- 
sociated with the resolution limit (0.2 AU). The disky shape 
is maintained by the high spin angular momentum of each 
clump. Depending on the mass and size, the total angular 
momenta of the clumps range from lO*'^ to 10^^ g cm'^/s. 
Assuming the total angular momentum is conserved during 
the evolution and a dopting a typical rad ius of a brown dwarf 
to be 7.0 X 10^ cm ([Burrows et al . 1998), the rotation speed 
at the edge of the brown dwarf would be 10 times the break 
up speed. Thus the proto-brown dwarfs disks must shed ro- 
tation in order to collapse down to a proto-BD. This result is 
consistent with recent mid-IR, sub-millimeter, and millime- 
ter observations of the spectral energy distributions (SEDs) 
of young brown dwarfs, which indicate that most of them 
are s urrounded by circumstellar disks (e.g., iMuench et al.l 
I2OOII ). A T Tauri accretion phase has also been indicated in 
some young B P disks by the detection o f broad, asymmet- 
rical Ha lines l|javawardhana et al.ll200^ ). The IR excess in 
SEDs was found to decrease as the age of BDs increases, 
which suggests that the BD disks evol ve from flared t o flat 
geometry similar to proto-stellar disks (|Mohantvll2005l ) . The 
clumps in our simulations are flat partly due to the assumed 
isothermal EOS. As explored in section [G] these dense con- 
densations are probably optically thick and would thus be 
hotter and more puffed up in nature. 

The simulations were terminated a few thousand years 
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Figure 9. A resultant clump with typical disky shape. Both pan- 
els are 14 AU on a side. The clump has disk radius about 4 AU. 
Left: The object projected in the x-y plane; Right: The object 
projected in x-z plane. The gray scale indicates the density in 
logarithmic space in code units Mq/AU~'^. 

after first object formed due to computer time limitations. 
In some objects, like the one shown in Figure |§J by the time 
the simulation was halted, the gas with lower specific an- 
gular momentum had started to condense into the center 
and a density gradient had also been established. Because 
the gravitational softening (0.2 AU) is comparable to the 
size of most clumps, our ability to resolve the structure and 
follow the evolution of the proto-BD disks is limited. How- 
ever, it can be expected in later evolution a central proto 
BD will form and the surrounding disk will accrete onto it 
on a viscous timescale. 

In the majority of cases where clumps formed from 
fragmentation of disks, the spin directions of the proto-BD 
disks are similar to their parent proto-stellar disks. In the 
case where the shock layer fragments, however, the resultant 
proto-BD disk may spin in a direction different from both 
disks. For the proto-BD disks that are still orbiting the star, 
further interaction with their parent disk will influence the 
final spin. 

5.3 Binary Brown Dwarfs and Multiple Systems 

Typically where fragmentation occurred more than one ob- 
ject was produced. About 15 percent of the objects are in 
BD-BD binary systems. This fraction is consistent with the 
observed binary fraction (aro und 10 ± 10 % at the rang e 
0.01 to 0.1 Mq) (cf. Figure S iHubber fc Whitworthll2005l ). 
In most cases, the binaries formed from secondary fragmen- 
tation of very large disky clumps. Further numerical stud- 
ies with higher resolution are necessary to investigate this 
question in detail because, although the Jeans mass was 
well resolved before fragmentation and most of the clumps 
are resolved by at least 1000 particles, the Jeans masses 
of the clump themselves decrease so quickly with increas- 
ing density that they are only marginally resolved in our 
simulations. Without resolving the Jeans mass, it has been 
suggested that fragmentation could be artificially enhanced 
(Bate fc Burkcrt 1997.). However, given the very high spin 
rate of the clump s, secondary fragmenta t ion is likely to oc- 
cur, according to lMatsumoto fc Hanawal (|2003l ). 

Most of the BD binary systems were still bound to the 
central stars when the simulations were terminated. For ex- 
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ample, in case "ConfS", the binary system orbits the star 
at ~ 140 AU, with the two components (mass 19 Mj and 
13 Mj) orbiting each other in a period of around 30 years 
and with separation about 3 AU. Each component has a 
highly flattened shape. In future evolution, it is expected 
that some of these binaries will merge to the star, leave the 
system, or be ejected in close interactions with other objects. 
The surviving systems are likely to remain hierarchical star- 
BD multiples. Thus the collision mechanism may explain 
the origin of observed star-BD multiple systems such as GL 
569, where a recently confirmed BD triple ( GL 569B) or- 
bits the primary, an M2 .5V star (GL 569 A) (|Gorlova et al.l 
l2003l : ISimon et al.ll2006l '). The system is about 100-125 Myrs 
old so no BD disks were observed. However, in younger BD 
multiples accretion signals have been detected (|Kraus et al.l 
l2006l ). indicating accretion disks surrounding components of 
the multiple. 

5.4 Orbits and Further Evolution of the Resultant 
Clumps 

Objects in our simulations typically formed at least 50 AU 
away from the closest star. Approximately 60 percent of the 
objects formed were not bound to the system. The percent- 
age is higher for coplanar encounters, in which the shock 
layer fragmentation can take place relatively far out in the 
gravitational potential of the stars. The unbound objects 
may leave the whole system to become free floating BDs 
(or BD binaries) or low mass stars. For the objects that are 
bound to the star, the orbital properties vary from case to 
case. About 5 percent of the orbits have very large semi- 
major axes (> 500 AU) but are still bound to the system. 
Since the outer proto-stellar disk was dispersing in most en- 
counters, sub-stellar objects in large orbits were unlikely to 
accrete additional gas to become low mass stars. Hence these 
objects may remain as wide-separation BD companions , 
which have been detected in observations (|Gizis et al.ll200ll ). 
Where more than one object orbits the star on nearby orbits 
dynamical interactions can take place r esulting in the pref- 
erent ial ejection of lower mass objects l|Reipurth fc Clarke! 
120011 ). The 11 objects that were ejected to significant dis- 
tances during the simulations were lower in mass. However, 
the majority of the objects were unbound and would ulti- 
mately become free-floating. 

For the clumps that stay in the dense region of the 
proto-stellar disks (within 100 AU), the future growth and 
evolution depend on the rate that the gas accretes onto the 
objects. Under the assumption of e flicient gas accretion, as in 
the competitive accretion scenario (|Reipurth fc Clarkell200ll : 
iBate et al.ll2003l ). these objects may eventually become low- 
mass stars over timescales much longer than those simulated 
here. 



6 COOLING APPROXIMATIONS 

6.1 Optically Thick Approximation 

The simulated disks produced dense regions as they evolved 
and reached the point where they were no longer uniformly 
optically thin. The disk material itself tended to remain op- 
tically thin (r < 1) and was thus capable of efficient cooling 



but the fragments within the disk reached high densities 
where the radiative cooling times became very long. Thus 
on the timescale of these simulations, the interiors of the 
fragments were effectively adiabatic. 

In earlier wor k on proto-planetary disk fragmentation 
iMaver et al.ll20()3 ) , the approximation of converting the en- 
tire simulation to adiabatic gas was employed to test the 
behaviour in this regime. This corresponds to the extreme 
where cooling is completely absent. A less extreme test with 
strong suppressed cooling is to smoothly interpolate from 
efficient black-body cooling when the gas is optically thin 
to no cooling in the optically thick regime. Optically thin 
black-body radiative losses are given by, 

L = AKaT'^M (4) 

where L is the luminosity, k is the opacity, T is the tem- 
perature and M is the total mass. The opacity was taken 
to be 1 cm^ g~^ at T = 60K with a power-law variation 
with temperature as T^ '^. This is consistent with the dust 
opacity valu es calculated by d'Alessio et al., as tabulated in 
iMeiial (|2004l '). over the range of 10-300 K. 

We then took a late stage of a representative run 
"Conf30" at 15000 years, after fragmentation had occurred, 
and measured the optical depth to infinity in 12 directions 
(corresponding the faces of a dodecahedron) for each parti- 
cle. This allowed us to estimate the averaged effective optical 
depth, Teff = -log^ e"^'"'dr2/47r. We found that the ef- 
fective optical depth correlates fairly well with gas density in 
the regime, p = lO-i"^ to 10"^^ g cm , where the material 
transitions from optically thin to (reff ~ 0.1) to optically 
thick (Tcff ~ 10) as t^b = iQ^-i+osep _ rpj^^ deviation 
of the fit is 0.16 dec in r. Assuming isotropic emission, we 
can then estimate the fraction of the emitted light (given 
by equation |4]) that escapes the disk entirely. We then as- 
sume that the rest of the emission is absorbed on-the-spot, 
an approximation commonly used for optically thick lines in 
models of emission nebulae, so that the net radiative loss, 
Lp, from one particle is Lp = 4KaTpmpe~^''^^^''''\ where rup, 
Tp and Pp are the particle mass, temperature and density, 
respectively. In the limit of infinitesimal particle masses, this 
expression is proportional to the contribution to the surface 
fiux from a mass of gas, dm = pdl dA = dr/n dA. In the case 
of an infinite plane-parallel slab, when an analytic expres- 
sion for Tcff is used and the flux contribution is integrated 
from the surface to high optical depth the total surface flux 
asymptotically approaches the value of aTp. 

The assumption for the main set of simulations was that 
the disk temperatures were set by heating from the central 
star offset by radiative losses so that the net temperatures 
were a function of the distance to the nearest star. To include 
the effect of heating by the star, the governing equation of 
the particle energy per unit mass, is, 

^ = 4,cae-""«('''')(r(r,tar)' - T^t) - -V.vj^ + A.hock.p, 
at p 

where r(rstar) is the temperature determined from an equi- 
librium between stellar heating and radiative losses at the 
radius Tatar from the nearest star, — -^V.v is the compres- 
sive heating term and Ashock is the shock heating term. 
Thus the model now includes compressive heating and shock 
heating. Outside the intense collision environment, where 
these process are unimportant in relation to stellar heat- 
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ing, the disk temperatures relax to the temperature as 
a function of radius assumed in the preceding sections. 
Overall, this similar to t he approximation employed by 
IStamatellos fc Whitworttil ([2002), where local cooling was 
based on the amount of radiation escaping from regions with 
similar physical conditions in spherical collapse models. In 
the approximation used here the radiative loss rates were 
calibrated from similar disk cases. 

Our approximation retains thermal energy indefinitely 
in the denser gas and is thus conservative with regard to 
ensuring that the fragmentation is not a result of favourable 
cooling assumptions. Detailed radiative transfer calculations 
should give results between this app roximati o n an d the 
isothermal EOS. It has been arg ued bv lRafikovl()2007l l that 
the optically thick parts of proto-stellar disks are likely to 
be convective, in which case our neglect of radiation as a 
heat transport mechanism may not have a large impact on 
the disk structure. However, the cooling of the optically thin 
surface layers are approximately correct in this model. Full 
radiative transfer has been implemented within the Gaso- 
line code (Rogers, Shen & Wadsley, in preparation) and will 
be applied to proto-stellar disk collisions as part of future 
work. 




ISO_clump. 



6.2 The Optically Thick Run 

We re-ran the simulation we used to calibrate the effective 
optical depth-density relationship using the new cooling ap- 
proximation. Figure [10] compares the disks at 14,000 years 
in the locally isothermal run (Panels that are labeled with 
"ISO_l" and "ISO_2") with the ones in the new run (Pan- 
els labeled with "THICK.l" and "THICK_2"). Using the 
isothermal equation of state, one of the disks fragmented 
and formed ten bound objects, including two binary sys- 
tems. However, when using the optical thick approximation, 
the cooling is suppressed and the disk does not fragment im- 
mediately. Instead, it is the tidal structure induced by the 
gravitational resonance that condensed and formed a clump 
at about 410 AU from the star. At later time steps, the gas 
in the tidal tail dissipated and the objects orbit the star at 
about 350 AU. The clump had a mass of 0.075 Mq, close 
to the hydrogen burning limit. 

Compared to the objects that formed in the correspond- 
ing isothermal run, the clump had a more spherical shape. 
Nevertheless, the angular momentum of the object was still 
very high, about 7.5 x 10^" g cm'^/s. By the end of the simu- 
lation we did not see flattening as the clump could not cool 
with this approximation. If no angular momentum was lost 
the final object would spin at 30 times the break up speed. 
Hence we expect the clump to collapse to a disk similar to 
the ones formed in the isothermal run and that jets and 
fragmentation could occur in nature as before. 

Thus with heavy cooling suppression, beyond that ex- 
pected in nature, we retain our main result, i.e. there was 
still gas condensation and the formation of clumps due to 
encounters. The number and the shape of the clumps at 
the end of simulation was strongly affected. Exploring the 
effect of differing thermal regulation on encounter-induced 
fragmentation will be a target for future work. 
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Figure 10. A Q = 1.6 disk encounter with isothermal EOS vs. 
our optical thick approximation. Left: The run with isothermal 
EOS (top and middle panels) and the typical shape of a clump 
(bottom panel); Right: The simulation that used optical thick 
approximation. The grey scale under each panel gives the density 
ranges in units of Mq/AU^. 



7 DISCUSSION AND SUMMARY 

In this work, the possibility of encounter-induced disk frag- 
mentation as a mechanism for forming sub-stellar objects 
was investigated. The significant increase in resolution over 
previous work and the use of a realistic initial disk model 
ensures that the initial fragmentation to form clumps was 
physically modeled within the constraints of our thermal 
regulation assumptions. These simulations have shown that 
encounters can trigger the formation of brown dwarfs as well 
as low mass stars. 

We found that strong instabilities can be triggered in 
stable disks with initial Qmin < 2.1, by lowering the local Q 
value with large gas inflows during an encounter. A relatively 
steep surface density proflle (E oc r~^^^) and a relative high 
ambient temperature (40 K) were used in our simulations. If 
disks have shallower density profiles, as was fo und in some 
observations of early disks l| Andre et al.|[200ol ). and hence 
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have large portions of their mass in the cooler, outer region, 
the rate of fragmentation could be even higher. 

By studying collisions with random initial orientations, 
it was also found that retrograde encounters (disks with ret- 
rograde spins relative to the orbital angular momentum) 
more effectively transport angular momentum, triggers gas 
inflows and ultimately produce more clumps on average than 
the prograde ones (Figure[6l). These results differs from those 
reported for simulations of star-disk encounters which indi- 
cated that prograde encou nters were more ef f icient in trig- 
gering disk instability (e.g. iBofiin et al ] |l998l : lLodato" et al.l 
l2007t l. The difference may relate to having an extended per- 
turber. 

Dense shock layers may occur during coplanar encoun- 
ters and result in objects as seen in previous work. How- 
ever, our results indicate that only slow coplanar encoun- 
ters (cases "Conff" and "Conf2"), where the shock layer 
had time to accrete substantial mass, can undergo frag- 
mentat ion. This differs from the results of IWatkins et al.l 
l|f998bl ). where a large fraction of their objects formed in 
the fragmentation of shock layers in both coplanar and non- 
coplanar cases. This discrepancy may be due to the use of 
more unstable initial disks, numerical problems related to 
lower resolution or the absence of vertical structure in the 
disk models (Appendix|X]| . ffowever it should be noted that 
coplanar encounters are very rare if the disk orientations are 
random. 

The velocity constraint found in our simulations implies 
the number of BDs produced in encounters is a strong func- 
tion of the velocity dispersion in the young star cluster. How- 
ever, the numerical trend as a function of cluster density can 
be complicated. On one hand, more massive and denser clus- 
ters usually have larger velocity dispersions and hence likely 
to produce fewer sub-stellar objects per encounter. On the 
other hand, high number density will increase the encounter 
rate which could potentially increase the BD numbers. 

The encounter mechanism is complementary to other 
mechanisms for producing BDs. Proto-stellar disks in close 
proximity are naturally produced when star cluster forma- 
tion is simulated from cloud collapse. Simulations of mas- 
sive clouds typically lack the resolution necessary to form 
large stable disks or t o follow subsequent encounters in de- 
tail (|Bate et al.ll2003l ). Thus this mechanism is compatible 
with ideas such as competitive accretion and turbulent frag- 
mentation. The relative importance of these processes for 
forming proto-stars may influence the typical separation and 
orientation of proto-stellar disks and thus affect the starting 
point for encounter driven object formation. These statis- 
tics also depend the larger environment, such as the size, 
the mass and the amount of turbulent motion in the star 
forming cloud. 

A dynamical encounter origin for some BDs might ex- 
plain observed differences in binary statistics between stars 
and BDs in different cluster environments. Some observa- 
tions have indicated that the ra tio of BDs to s t ars is higher 
in denser clusters. For example. ISlesnick et al.l (j2004l ) found 
in Trapezium, the number ratio (Ras) of sub-stellar objects 
(0.02 £ M/Mq ^ 0.08) over sta r s (0.0 8 < M/Mq < 10) is 
about 0.20, while iLuhman et al] l|2003h calculated for Tau- 
rus aggregates and IC 348 Ras is about 0.f4 and 0.12 respec- 
tively, significantly lower than Trapezium. The results are 
not conclusive as other surveys have indicated that Ras in 



Taur us is comparable with that in Trapezium (jLevine et al.l 
|2006| . and references therein). To make further predictions 
about how the BD ratio may be related to the cluster proper- 
ties, more comprehensive simulations and parameter studies 
would be required, in which the encounter rate in different 
cluster environment is estimated in detail and the encounter 
velocity parameter space is well sampled. 

Almost all the clumps formed in the isothermal simula- 
tions have flattened shapes. In the future evolution, we ex- 
pect these clumps to centrally condense to form a proto-BD 
with a disk. Material should accrete onto the central object 
on the viscous timescale. Observations indicate that the ac- 
cretion rate for low mass stars and sub-stellar objects can 
be approximately fitted by M ~ 10"*Mq yr~^ (M/Mq)^ 
(MuzcroUo ct al. 200j). Hence for a relatively massive brown 
dwarf with 0.05 Mq the accretion rate is around 2.5 x 
10~^^ A/q yr~^. With this rate the BD disk can have 
a mass as low as 10"'* A/q and still have lifetime com- 
parable to proto-stellar disks. To shed angular momen- 
tum, a proto-BD disk might also launch bipolar outflows or 
form sub-stellar companions. Since most BDs formed dur- 
ing encounter-induced disk fragmentation have large excess 
angular momenta, we expect outflows or companions to be 
be common for young BD formed this way. Recent obser- 
vations have found forbi dden emission lin es in BD spec- 
tra, suggesting outflows l|Natta et al.ll2004l ). In particular, 
an outflow fro m a young BDs p Oph 102 has been spa- 
tially resolved (jWhelan et al.|[2005l ). Some of proto-BD disks 
in our simulations underwent secondary fragmentation and 
form binary BD systems (cf . Section 15. 3p . Due to resolu- 
tion and physical limitations, this result needs to be con- 
firmed with better simulations. How ever, the rapid rotation 
of th ese disk favours fragmentation ()Matsumoto fc Hanawal 
120031 ) which can form companions of a range of masses. This 
provides a mechanism for forming the observed BD-BD bi- 
naries with the small separations pea king at ^ 1 to 4 AU and 
with almost equal mass components (IWhitworth et al]|2007l . 
and references therein). Planetary-mass companions have 
also been observed surrounding young BD (e.g. 2M1207b 
IChauvin et "al]|2005l ). 

Possible future work could follow several directions. 
Firstly, the predictive power would be improved if parame- 
ters other than disk configuration were to be well sampled, 
particularly the velocity of the encounter, the disk properties 
and the star masses. Combined with a detailed calculation 
of encounter rates as a function of cluster density and ve- 
locity dispersion, quantitative results could be obtained to 
directly compare with the observations of BD to star ratios, 
Ras- Secondly, when the optical-thin approximation fails in 
high density regions, the effects of shock heating and ra- 
diative cooling become more important and should be fur- 
ther investigated. Including radiative transfer in the simu- 
lation is likely to reduce the number of objects formed per 
encounter, and provide for more realistic evolution of the 
proto-BD disks (as indicated in the optically thick approxi- 
mation, section IS}. Alternately, with simplified models, the 
evolution of the clump themselves could be followed over 
long timescales. 



16 S. Shen et al. 



APPENDIX A: RESOLUTION 
REQUIREMENTS 

When designing the sim ulation the resolu t ion r equirement 
was estimated from the iBate fc Burkertl criterion, 
i.e., the Jeans mass should be resolved by at least twice the 
number of neighbour particles, N Neighbor - In our simulations 
we used exactly 32 neigbours. In the initial condition, the 
density of our disk varies from 10""5/cm^ to 10"^°g/cm^. 
Using our temperature profile, described in section [J] the 
Jeans mass at the highest density for the initial disk was 
estimated to be 0.0095 Mq. Hence resolving the initial disk 
of 0.6 Mq, which contains about 63 of these smallest Jeans 
masses requires about 4000 particles for each disk. During 
the encounter, however, as shown in Section |31 large inflows 
of gas and unstable modes increase the density significantly. 
The density typically increases by three orders of magnitude 
locally where fragmentation will then occur. Therefore to re- 
solve the smallest Jeans mass during the simulation signifi- 
cantly more particles were necessary. To correctly simulate 
gas with densities up to 10~^°(?/cm^ about 127000 particles 
were needed. We used 200000 particles for each disk, which 
met the above requirement and ensured that the fragmen- 
tation was physically modeled. 

The effects of resolution and the finite scale height were 
studied with two-dimensional simulations and through an- 
alytical arguments by N elson C200(|). We performed two 
sets of comparison simulations to investigate these effects 
in three dimensions. The first set consisted of three simu- 
lations of isolated disks with the same surface density and 
temperature profile and a minimum Toomre Q of 1.3, well 
above the critical value of 0.8. The first disk was modeled 
by the single layer of 2000 particles. The second had 20000 
particles but still had no vertical structure. The third one 
used 200000 particles and the vertical structure was mod- 
eled using the methods described in section (2] Figure lAll 
depicts the results for each disk at a later time step. 

While no structure formed in the third, well-resolved 
simulation, strong gravitational instabilities and fragmenta- 
tion occurred in the first two within several outer orbital 
periods. Note that in the first case with 2000 particles the 
Jeans mass was not resolved in the high density region of the 
initial disk. For the second case although the Jeans mass was 
fully resolved everywhere in the initial condition (not nec- 
essary during the simulation) , the lack of vertical structure 
in the model was still destabilizing and artificial fragmenta- 
tion occurred. Although it is not clear from these two simu- 
lations which factor (resolving Jeans mass or modeling the 
finite thickness) is more important to prevent unphysical 
fragmentation, the test simulations show that it is essential 
to consider both factors in the simulation. 

The second series consist of two simulations of retro- 
grade dynamical encounters of disks with Q = 2.1. The disks 
were very stable in isolated runs. The dynamical parameters 
of the encounter was the same as case "Q3" in Table [T] In 
the first simulation the disks were modeled by single layer 
of 2000 particles, while in the second the disk was resolved 
by 200000 particles and had vertical structures. The results 
is depicted in Figure 1^2] 

Similar to the isolated runs, without properly resolv- 
ing the Jeans mass and scale height, the encounter triggered 
fragmentation at time step of 19700 years in the shock layer 



between two coplanar disks (Panel 2, Figure R2)) . In the sec- 
ond simulation, although the shock layer was still produced 
(Panel 4), no fragmentation occurred in the layer. The en- 
counter was also modeled using 400000 particles with the 
same result. 

In summary, regardless of whether the disks are pas- 
sively evolving or undergo a physical impact, numerical ef- 
fects may be significant if the Jeans mass is not resolved or 
if no vertical structure is presented in the disk model. The 
major numerical effect is then enhanced instability and the 
formation of artificial clumps. Our resolution of 200000 par- 
ticles was chosen based on a conservative estimate of resolu- 
tion requirements, e.g., it assumed the density of the whole 
disk could increase up to 10~^^(;/cm^. In fact, most outer 
region of the disk had density below 1Q~^^ g / cm?' when frag- 
mentation initial occurred. With this resolution it was en- 
sured that the initial disk fragmentation (if any) was phys- 
ical. We also performed simulations at higher resolution 
(400000 and 800000 particles) for some specific cases and 
the results converged. 
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Figure Al. Isolated disk simulations with different resolutions and disk models. All three disks have Q = 1.3 initially. Time steps for 
each snapshot are given in the upper left corners. Boxed regions indicate the fragments (if any). The grey scale under each panel gives 
the density ranges in units of Mq/AU^ . 
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Figure A2. Disk encounter simulations with different resolutions and disk models. The disks have Q = 2.1 initially. Time steps for each 
snapshot are given in the upper left corners. The first two panels (panel 1-2) depicted the low resolution run and panel 3-4 are from the 
higher resolution one. Boxed regions indicate the fragments (if any). The grey scale under each panel gives the density ranges in units of 
Mq/AU^. 
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